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Abstract 

OO ' 

The relativistic four-quark equations in the framework of the dispersion relation tech- 
nique are derived. The four-quark amplitudes of the scalar tetraquarks with open 
bottom, including u, d, s and bottom quarks, are constructed. The poles of these 
amplitudes determine the masses and widths of scalar tetraquarks. 
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I. Introduction. 

The observation of the X(3872) [1-4], the first of the XYZ particles to be seen, 
brought forward the hope that a multiquark state has been received. Belle Collabora- 
tion discovered the X(3940) in the reaction e + e~ — > J/tp + X [5]. The fact that the 
newly found states do not fit quark model calculations [6] has triggered the interest 
in the charmonium-like and the charmed states. Maiani et al. advocate a tetraquark 
explanation for the X(3872) [7, 8]. Ebert et al. [9] calculated that the light scalar 
tetraquark lies above the open charm threshold and is broad, while Maiani et al. ob- 
! tained [10] that this state lies a few MeV below this threshold. In our paper [11] 

the tetraquark X(3700) with the spin-parity J pc = ++ decays to rjrj c . The calculated 
width of this state is equal to r ++ = 140 MeV. 

In the present paper the relativistic four-quark equations are found in the framework 
of coupled- channel formalism. The dynamical mixing between the meson-meson states 
and the four-quark states is considered [12 - 14]. The masses and the widths of the 
scalar tetraquarks with open bottom are calculated. 

II. Four-Quark Amplitudes for the Tetraquarks with Open Bottom. 

We derive the relativistic four-quark equations in the framework of the dispersion 
relations technique. 

The correct equations for the amplitude are obtained by taking into account all 
possible subamplitudes. It corresponds to the division of complete system into subsys- 
tems with the smaller number of particles. Then one should represent a four-particle 
amplitude as a sum of six subamplitudes: 

A = A 12 + A 13 + A u + A 23 + A 2i + A u . (1) 

This defines the division of the diagrams into groups according to the certain pair 
interaction of particles. The total amplitude can be represented graphically as a sum 
of diagrams. 
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We need to consider only one group of diagrams and the amplitude corresponding, 
for example A 12 . 

The relativistic generalization of the Faddeev-Yakubovsky approach [15, 16] for the 
tetraquark is obtained. We shall construct the four-quark amplitude of buuu tetraquark 
with the spin-parity J pc = ++ in which the quark amplitudes with quantum num- 
bers of 0~ + and 1 mesons are included. The set of diagrams associated with the 
amplitude A\ 2 can further broken down into four groups corresponding to subampli- 
tudes: Ai(s,s 12 ,su), A 2 {s, s 23 , %), A 3 (s, s 2 3, s 123 ), A 4 (s, s 14 , s 124 ), if we consider the 
tetraquark with the spin-parity J pc = ++ (buuu). 

Here s ik is the two-particle subenergy squared, corresponds to the energy 
squared of particles i, j, k and s is the system total energy squared. 

In order to represent the subamplitudes Ai(s, Si 2 , s 34 ), A 2 (s, s 23 , s 14 ), A 3 (s, s 23 , S123) 
and A 4 (s, Si 4 , si 24 ) in the form of dispersion relations it is necessary to define the ampli- 
tudes of quark-antiquark interaction a n (si k ). The pair quarks amplitudes qq — > qq are 
calculated in the framework of the dispersion N/ D method with the input four-fermion 
interaction [17 - 19] with quantum numbers of the gluon [20]. The regularization of 
the dispersion integral for the D-function is carried out with the cutoff parameter A. 
The four-quark interaction is considered as an input [20]: 

9v {q^Ifl^q) 2 + 9 ( v (q~Mn»q) (sXl^s) + g ( y s) (sX-y^s) 2 . (2) 
Here // is the unity matrix in the flavor space (u,d). A are the color Gell-Mann 
matrices. Dimensional constants of the four-fermion interaction gv, g$ and gy^ are 
parameters of the model. At gy = gy^ = gy^ the flavor SU (3) j symmetry occurs. The 
strange quark violates the flavor SU(3)f symmetry. In order to avoid an additional 
violation parameters, we introduce the scale shift of the dimensional parameters 



m 2 (m + m s ) 2 {s) _ m 2 {ss) 
9--I9V- ^5 9v -~l9v ■ (3) 

A = ^) (4) 
(rrii + m k y 

Here m« and are the quark masses in the intermediate state of the quark loop. 
Dimensionless parameters g and A are supposed to be constants which are independent 
of the quark interaction type. The applicability of Eq. (2) is verified by the success of 
De Rujula-Georgi-Glashow quark model [21], where only the short-range part of Breit 
potential connected with the gluon exchange is responsible for the mass splitting in 
hadron multiplets. 

We use the results of our relativistic quark model [20] and write down the pair 
quarks amplitude in the form: 

a n {sik) = 5-7 — r , (5) 

(m i +m fc ) 2 A 

B n (s lk )= f ^ P " (4)G " (4) ■ (6) 

, J , a S' lk -S lk 

Here G n (s ik ) are the quark-antiquark vertex functions. The vertex functions are de- 
termined by the contribution of the crossing channels. The vertex functions satisfy 
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the Fierz relations. All of these vertex functions are generated from gy, 9v^ an d 
gy S \ B n (sik), p n ( s ik) are the Chew-Mandelstam functions with cutoff A and the phase 
spaces, respectively. 

Here n — 1 determines a gg-pairs with J pc = 0~ + in the l c color state, n = 2 
corresponds to a gg-pairs with J pc = 1 in the l c color state, and n = 3 defines the 
gg-pairs corresponding to tetraquarks with quantum numbers J pc = ++ . 

In the case in question, the interacting quarks do not produce a bound state; there- 
fore, the integration in Eqs. (7) - (10) is carried out from the threshold (m 8 + m^) 2 
to the cutoff A(ik). The coupled integral equation systems (the tetraquark state with 
n = 3 and J pc = ++ buuu) can be described as: 



^iB 2 (s 12 )B 2 (s 34 ) 
[1 - B 2 (s 12 )][l - B 2 (s u )} 



Ai(s, S12, s 34 ) — — -!— ^ — 12 2 ^ 4 — — + 2J 2 (s 12 , S34, 2, 2)A 3 (s, s 23 , s' 123 ) 



+ 2J 2 (s 12 ,S34,2,2)A4(s,s' 14 ,s' 124 ), (7) 

/I f ^ A 2j Bi(s 23 )Si(Si 4 ) 1 i\4 / ' / \ 

A-2{S, S 23 , S 14 ) - — Bi ( S23 ^i _ B^Jsu)] 2l - S23 ' Sl4 ' ' L ) A ^ S ' S 34? S 234J 

+ 2J 2 (s 2 3,Si4,l,l)A 4 (s,s' 12 ,s' 12 3) , (8) 

A3 ^3 ( 5 23 ) A it ^ 

M{S, S 23 , S123) = Tj 5-7 77 + 2 J 3 (S 23 , 3)Ai(s, S 12 , S34) + J 3 (S23, 3)A 2 (s, s 12 , S34) 

+ ^1(^23, 3)A 4 (s, s 34 , s 234 ) + ^i(s23 5 3)A 3 (s, s' 12 , s' 123 

), (9) 

A 4 (s, su, S124) = Fr^%~rSi + 2 ^s(si4, 3)Ai(s, s' 13 , s' 24 ) + 2J 3 (s 14 , 3)A 2 (s, s' 13 , s 24 ) 

[1 - # 3 (Si 4 JJ 

+ 2J 1 (s 14 , 3)A 3 (s, s' 14 , s' 134 ) + 2J 1 (s 14 , 3)A 4 (s, 

S 14> S 134) ) 

(10) 

where Aj, i = 1,2,3,4 are the current constants. They do not affect the mass spectrum 
of tetraquarks. We introduce the integral operators: 

(m 1 +m 2 ) 2 A ^ 

f(o n - G i( s ^) / ds' 12 Gi(s' 12 )pi(s' 12 ) f dz x . . 

JlKSi2,l) - 7- —, 77 / ; / — , [ii 

[l-Bi(s 12 )\ / J it s' 12 -s 12 J 2 

(ml+m.2) — 1 



^2(Sl2,S 3 4,/,P) = 



Gi(s 12 )G p (s u ) f ds' 12 Gi(s' 12 )pi(s' 12 ) 



[1 -Bi(s 12 )][l -B p (s 34 )] no ^* s' 12 -si2 

(mi+m2j / 



(m3 + m 4 ) 2 A 



/" rfs 34 G'p(s3 4 )p p (s3 4 ) Vdz 3 + r 1 dz 4 , . 

J tt sL-S3 4 J 2 7 2 ' 



(m3+ni4) 2 34 -1 
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J , « gi^f^A) J_ J ds' 12 Gi(s' 12 ,A)pi(s' 12 ) 

3lSl2 ' J ~ [1- 5,(3X2, A)] 47T 7 7T '~ S' 12 -S 12 '~ 

(mi+m 2 ) i 

+1 , +i 4 

x /" ^ / ' dzjdz 2 . (13) 



-l -l 



J I — z 2 — z\ — z\ + 2zz\ z 2 



here /, p are equal to 1 — 3. 

In Eqs. (11) and (13) z\ is the cosine of the angle between the relative momentum 
of particles 1 and 2 in the intermediate state and the momentum of the particle 3 in 
the final state, taken in the cm. of particles 1 and 2. In Eq. (13) z is the cosine of the 
angle between the momenta of particles 3 and 4 in the final state, taken in the cm. 
of particles 1 and 2. z 2 is the cosine of the angle between the relative momentum of 
particles 1 and 2 in the intermediate state and the momentum of the particle 4 in the 
final state, is taken in the cm. of particles 1 and 2. In Eq. (12): z 3 is the cosine of the 
angle between the relative momentum of particles 1 and 2 in the intermediate state 
and the relative momentum of particles 3 and 4 in the intermediate state, taken in the 
cm. of particles 1 and 2. z 4 is the cosine of the angle between the relative momentum 
of particles 3 and 4 in the intermediate state and that of the momentum of the particle 
1 in the intermediate state, taken in the cm. of particles 3 and 4. 

We can pass from the integration over the cosines of the angles to the integration 
over the subenergies [22]. 

Let us extract two-particle singularities in the amplitudes Ai(s, s± 2 , S34), A 2 (s, s 2 3, S14), 
A 3 (s, S23, S123) and A 4 (s, s u , s 124 ): 

Af x _ ai(s,s ik ,s lm )B 2 (s ik )B 2 (s lm ) 

A(s ' s * Slm) ~ FW^« ' ( } 

. ( x _ o^js, Si k) si m )Bi(s ik )Bi(si m ) 

A 2 \S, S ik , Si m ) — — r, 5-7 xTj-j 7T7 ^~ 1 

[1 - 5i(s ifc )J[l - Bi{si m )\ 

A j (s,s ik ,Siki) = 1 _ B ^ s .^ > J = 3-4. (16) 

We do not extract three-particles singularities, because they are weaker than two- 
particle singularities. 

We used the classification of singularities, which was proposed in paper [23]. The 
construction of the approximate solution of Eqs. (7) - (10) is based on the extraction of 
the leading singularities of the amplitudes. The main singularities in s ik ~ (m^ + m k ) 2 
are from pair rescattering of the particles % and k. First of all there are threshold square- 
root singularities. Also possible are pole singularities which correspond to the bound 
states. The amplitudes apart from two-particle singularities have triangular singulari- 
ties and the singularities defining the interactions of four particles. Such classification 
allows us to search the corresponding solution of Eqs. (7) - (10) by taking into account 
some definite number of leading singularities and neglecting all the weaker ones. We 
consider the approximation which defines two-particle, triangle and four-particle singu- 
larities. The functions a^s, s 12 , s 34 ), a 2 (s, s 23 , s 14 ), a 3 (s, s 23 , s 123 ) and a 4 (s, s u , s 124 ) 
are the smooth functions of s ik , s ik i, s as compared with the singular part of the ampli- 
tude, hence they can be expanded in a series in the singulary point and only the first 
term of this series should be employed further. Using this classification, one defines 
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the reduced amplitudes a., a 2 , a 3 , a 4 as well as the £?-functions in the middle point 
of physical region of Dalitz-plot at the point s : 



4 fe = 0.25(m; + m k ) 2 s , 
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3 3 s + 2 £ m\ 

s 123 = 0.25s (mi + m k ) 2 -J2 m2 i , s = i — — . (17) 

%f i=1 0.25 E (mi + m fc ) 2 

Such a choice of point so allows us to replace integral Eqs. (7) - (10) by the algebraic 
equations (18) - (21) respectively: 

ai = Ai + 2a 3 J5!(2, 2, 3) + 2a A JB 2 (2, 2, 3) , (18) 
a 2 = \ 2 + 2a 3 JS 3 (1, 1, 3) + 2a 4 J5 4 (l, 1,3), (19) 
a 3 = A 3 + 2ai JCi(3, 2, 2) + 2a 2 JC 2 (3, 1, 1) + a 4 ^i(3) + a 3 JA 2 (3) , (20) 

a 4 = A 4 + 2ai JC 3 (3, 2, 2) + 2a 2 JC 4 (3, 1, 1) + a 3 JA 3 (3) + a 4 J7L 4 (3) . (21) 

We use the functions JAi(l), JBi(l,p,r), JCi(l,p,r) (l,p,r = 1 — 3), which are 
determined by the various Sg fc (Eq. 17). These functions are similar to the functions: 

j a m Gf(sl 2 )Bf(sf) f ds' 12 Pl (s' 12 ) + r 1 dz 1 1 
M ) ~ B l (s?) J * s> 2 -s 12 J 2 ( j 

(mi+m2)' ! —J- 



JB.il, p,r) = 



Bi(Sq 2 )B p (Sq 4 ) J 7T S' 12 -S 12 

(mi+m2)^ 



(m3+m4) 2 A 

'-1 +1 



x 



? rf^4 Pp(4J /• dz 3 f dz 4 1 

J 7T sL-smJ 2 J 2 l-BJs'n,)' 1 ' 



■S 34 — S 34 J I J I 1 — -D r 1^23 J 
(m 3 +m4j' ! -1 



(j7l[ +^2)^^ 



JC 3 (l,p,r) = 



GUsl 2 ,A)B p (sf)B r (s^)l-B l (sl 2 ) 1 > ds' 12 Pl {s' 12 ) 



1-^(4 2 ,A) ^(SO 2 ) 47T 7 7T s' 12 -s 12 

(mi+m2]' : 



+1^ +1 4 1 
x / — — I dz I dz 2 - 



1 2 _ 1 - ,2 2 - ^ - z§ + 2^2:^2 

Z 2 V 



x [l-s p ( s23 )][l-s r ( s ' 14 )]' (24) 
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A(ik) 



A(ik), 



if A(^)<(y% + m 3 f 



(25) 



(v^sim + m 3 ) 2 , 



if A(ik) > (^ + m 3 ) 2 



The other choices of point s do not change essentially the contributions of a\, a 2 , 
a 3 and a 4 , therefore we omit the indices s 1 q. Since the vertex functions depend only 
slightly on energy it is possible to treat them as constants in our approximation. 

The solutions of the system of equations are considered as: 



where zeros of D(s) determinants define the masses of bound states of tetraquarks. 
Fi(s, Aj) determine the contributions of subamplitudes to the tetraquark amplitude. 

III. Calculation results. 

Our calculations do not include the new parameters. We use the cutoff A = 7.63 and 
the gluon coupling constant g = 1.53, which are determined by fixing the tetraquark 
masses for the states with the hidden bottom [24, 25]. The widths of the tetraquarks 
with the open bottom are fitted by the fixing width T 2 ++ = (39 ± 26) MeV [26] for 
the S'-wave tetraquark with the hidden charm X(3940). The quark masses of model 
m u,d — 385 MeV, m s = 510 MeV and mj = 4787 Me V coincide with our paper ones 
[25]. The masses and the widths of meson- meson states with the spin-parity J pc = ++ 
are given in Table I. In our paper we predicted the tetraquark (buuu) with the mass 
M = 5914 MeV and the width r ++ = 104 MeV. We calculated the masses of X 6 (6020) 
M = 6017 Me V and the width r ++ = 69 MeV (channels B Q s r] and B + K~). The 
tetraquark (bsus) have the mass M = 6122 MeV and the width T ++ = 48 MeV. The 
tetraquarks with the open bottom and the spin-parity J pc = 1 ++ , 2 ++ have only the 
weak decays. 

The functions Fj(s,Aj) (Eq. (26)) allow us to obtain the overlap factors / for 
the tetraquarks. We calculated the overlap factors / and the phase spaces p for the 
reactions X —> M X M 2 (Table I). We considered the formula r ~ f 2 x p [27], there p is 
the phase space. 

In the open bottom sector the scalar tetraquarks have relatively small width ~ 
50 — 100 MeV, so in principle, these exotic states could be observed. 
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a i (s)=F i (s,X i )/D(s), 



(26) 
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Table I. Masses, widths, overlap factors / and phase spaces p of scalar tetraquarks with 
open bottom. 



Tetraquark 


(channels) 


£ 

J 


P 


Mass {MeV) 


Widths {MeV) 


X fe (5910) 


B + T] 


0.54 


0.103 


5914 


104 


X 6 (6020) 


B+K- 


0.32 
0.23 


0.108 
0.171 


6017 


69 


X fe (6120) 


B»K+ 


0.283 


0.174 


6122 


48 
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